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A circular magnetic disk with biaxial magnetocrystalline anisotropy has four stable magnetization states 
which can be used to encode a pixel's shade in a black/gray /white image. By solving the Landau-Lifshitz- 
Gilbert equation, we show that if moderate noise deflects the magnetization slightly from a stable state, it 
always returns to the original state, thereby automatically de- noising the corrupted image. The same system 
can compare a noisy input image with a stored image and make a matching decision using magneto-tunneling 
junctions. These tasks are executed at ultrahigh speeds (~2 ns for a 512x512 pixel image). 



Consider a single crystal single-domain nanomagnet in 
the shape of a circular disk that has no in-plane shape 
anisotropy but has biaxial magnetocrystalline anisotropy. 
This magnet has four stable magnetization directions 
('up', 'right', 'down', 'left') as shown in Fig. 1 (a)^. 

An array of such magnets can store a black-gray- white 
image by storing the shade of every pixel in one of the 
four stable states of a magnet. If corrupted by moderate 
noise, the magnetization dynamics automatically recover 
the stored image, so that there is built-in error correction. 
When integrated with magneto-tunneling junctions, the 
same array can recognize images by comparing them with 
stored images pixel by pixel. These tasks are executed 
with ultrahigh speed since no software is needed. 
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FIG. 1. (a) A circular disk-shaped single-crystal Ni layer with 
biaxial anisotropy, creating four possible magnetization di- 
rections or 'easy' axes - 'up', 'down', 'left' and 'right'. The 
energy landscape is depicted by the saddle-shaped curve, re- 
sulting in the four energy minima along the axes, (b) The 
spherical coordinate system used to define the magnetization 
vector and the color scheme used to designate quadrants, (c) 
Color plot of the magnetization vector's final state for all al- 
lowed values of the initial state {0 initial A initial)- 
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The four states of the magnet encode black, gray, white 
and gray. We could have encoded two different shades of 
gray in the four states, but the 3-shade scheme is sim- 
pler. Noise can corrupt the stored image by perturbing 
the magnetization vector and deflecting it away from the 
initial stable state. The ensuing magnetization dynam- 
ics is studied using the Landau-Lifshitz-Gilbert (LLG) 
equation to determine the final state of every perturbed 
magnet (or pixel). For this purpose, we assume that the 
magnets are spaced far enough apart that dipole inter- 
action between them can be ignored. Since the useful 
pixel density in any optical image is limited by the wave- 
length of light (0.5 - 1.0 jiia for visible light), we can 
space the magnets apart by 1 /im, which makes the inter- 
magnet dipole interaction negligible. Using knowledge of 
the magnetization dynamics of such a 4- state magnet, we 
demonstrate image recovery in a 512x512 pixel black- 
white-gray image in ~2 ns while conventional field pro- 
grammable gate array (FPGA)- and CMOS-based image 
filters would have taken several ji^. 

In order to demonstrate image recovery, we consider 
an isolated circular nanomagnet of nickel with diameter 
100 nm and thickness 10 nm. Nickel has an 'easy' axis 
of magnetization along the (111) direction, a 'medium' 
axis along the (110) direction and a 'hard' axis along the 
(100) direction. We assume the nickel layer to be in the 
(001) plane. In this plane, the 'easy' axis is the (110) di- 
rection ('medium' for the entire crystal) while the 'hard' 
axis is along the (100) direction. Because the thickness 
of the magnet is ten times smaller than the diameter, 
out-of-plane excursion of the magnetization vector is en- 
ergetically costly, albeit not impossible. 

As illustrated in Fig. 1 (a), the easy axes of a single- 
crystal Ni nanomagnet in the x-y plane (magnet's plane) 
are along the [110], [1 10], [110] and [110] directions, in 
Miller notation. The saddle-shaped curve and arrows de- 
pict the potential energy landscape and stable directions 
of the magnetization vectors in the four (degenerate) 
minimum energy locations. The coordinate axis system 
is rotated by 45° in the magnet's plane (i.e. (j) ^— 7r/4), 
where (j) is the angle subtended by the in-plane compo- 
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nent of the magnetization vector with the +x-axis. This 
rotational transformation ensures that the energy min- 
ima are along the ±x-axes ((/) = 0°or±180°) and ±y-axes 

= ±90°). The angle subtended by the magnetization 
vector with the z-axis is 6. When 6 = 90°, the magneti- 
zation vector is in the plane of the magnet. 

If noise deflects the magnetization vector of a mag- 
net away from its initial stable state to a new state, the 
latter evolves with time in accordance with the Landau- 
Lifshitz-Gilbert (LLC) equation. We take the new state 
as the initial condition and solve the LLG equation to de- 
termine the final state. We find that as long as the noise 
amplitude is not large enough to defiect the magnetiza- 
tion vector closer to another stable orientation or cause 
large out-of-plane excursion, the vector always returns to 
the initial stable orientation. This happens relatively fast 
in a few ns (for realistic parameters). Thus, the system 
acts like associative memory and could be used in that 
role as well. However, here we study image recovery. 

The total magnetic energy of the single-domain nano- 
magnet per unit voulme is 
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where Ki and K2 are the first and second order magne- 
tocrystalline anisotropy constants, /io is the permeability 
of free space, Ms is the saturation magnetization and 
Nx^Ny^Nz are the components of the demagnetization 
factors along the x-, y-, and z-axes, respectively. Since we 
consider a circular nanomagnet, = Ny. These demag- 
netizing factors are determined from the oblate spheroid 
estimation^. 

Since the magnet is single-domain, the magnetization 
M has a constant magnitude. This allows us to define 
a unit vector 11^ = M/|M| = in the radial direction. 
The magnetization at any instant is therefore uniquely 
specified by 0{t) and (j){t). 

In any arbitrary orientation [0{t)^ 4^(f)\i the magnetiza- 
tion vector experiences a torque per unit volume due to 
magnetocrystalline and shape anisotropy. The latter can 
be written as 

T^(t) = -n^ X VEtotalit) = -er X VEtotal {0{t),^{t)) , 

where Etotai (t) is given by Equation M . It can be verified 
that the torque vanishes when = 11/2 and (j) = nir /2 [n 
= 0, 1, 2, 3], indicating that these are the four stable ori- 
entations. Whenever the magnetization is defiected from 
a stable orientation by noise, it experiences the above 
restoring torque that attempts to return the magnetiza- 
tion to the initial state. 

In order to study the restoring dynamics, we solve the 
Landau-Lifshitz-Gilbert (LLG) equation 
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where a is the Gilbert damping coefficient (0.045 for Ni), 
Q is magnet's volume, and fiB = Bohr magneton. 

From the last three equations, we derive two coupled 
equations that describe the time-evolution of 0{t) and 
<p{t) (see the accompanying supplementary materiaP). 
After noise deffects the magnetization to a new orien- 
tation, we take that new orientation as the initial ori- 
entation [Oinitiah(t^initiai] and solvc the above coupled 
equations numerically to determine the final state of the 
magnetization vector [Ofinah^* final]- The final state is 
reached when the orientation does not change further. 

In our simulations, the initial conditions span the space 
of (/^initial from 0° to 360° (with four narrow excluded 
regions; see later) and OinUiai from 85° to 95°. The an- 
gular resolution is 1°. The results outside the interval 
[85° < Oinitiai ^ 95°] are not considered since they corre- 
spond to large out-of-plane excursions of the magnetiza- 
tion vector that can only occur when the noise amplitude 
is extremely high. Such large perturbations induce non- 
linear effects and will be discussed elsewhere. 

The narrow regions excluded from the (j)initiai space 
are the following. The states (j) = 45°, 135°, 225° and 
315° are the highest energy states when = 90°. If the 
magnetization is driven to any of these states by noise 
(unlikely since these are the highest energy states), it 
has equal probability of decaying to either one of its two 
neighboring minimum energy (stable) states. Therefore, 
we exclude a region of 1° around these critical points. 

In order to convey the results of magnetization dynam- 
ics obtained by solving the LLG equation, we use the fol- 
lowing color scheme (see Fig. 1(b)): blue corresponds to 
(/) = 0°, red to (j) = 90°, yellow to (j) = 180°, and green 
to (j) = 270°. Fig. 1(c) illustrates the magnetization 
vector's final orientation for all allowed values of initial 
(perturbed) orientations Oinuiai and (t)initiai- If Oinitiai 
and (^initial are in the blue region, then the final state is 
always (j) final = 0° and final = 90°. Similarly, if OiniUai 
and (^initial are in the red region, then the final state is al- 
ways (j) final = 90° and 6 final = 90°, etc. This plot shows 
that if Oinitiai = 90°, then the final state for the range 
of initial states -44° < ^initial < 44° is (j) final = 0° and 
final = 90°. Similarly, the final state for Oinitiai = 90° 
and the range of initial states 46° < (/^initial < 134° is 
4^ final — 90° and final = 90°, and so on. This means 
that if noise defiects the magnetization vector from a sta- 
ble state while still staying in-plane, the vector always re- 
turns to the original state as long as the angular in-plane 
defiection is within ±44°. 

However, when Oinitiai = 89°, i.e. noise causes -1° 
out-of-plane excursion of the magnetization vector, the 
regions shift by -5°. Now, if the magnetization's ini- 
tial state was in the interval [-50° < (j^initiai 

< 40°], the 



the final state will be 



0°, and so on. This 
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5° shift happens because the out-of-plane excursion of 
the magnetization vector generates an additional pre- 
cessional torque acting on the magnetization vector ow- 
ing to the coupled — (j) dynamics. Similarly, when 
Oinitiai = 91°, a shift of +5° takes place. As long as 
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FIG. 2. (a) A pixel's shade is given a value of (black), 0.5 
(gray) or 1 (white), (b) Original 512x512 pixel image. Each 
pixel's shade is encoded into the magnetization state of the 
corresponding magnet, (c) The eye of the bird is blown up for 
clarity, (d) Noise corrupts the image at time t = by deflect- 
ing the magnetization vectors from their initial states. The 
out-of-plane excursion is restricted to 1°. (e) The partially 
recovered image at t = 1 ns. (f) The fully recovered image at 
t = 2 ns. 



the out-of-plane excursion AOinitiai{= initial - 90°) is 
small, the shift is linearly proportional to AOiniUai- For 
^^initiai > 5°, the shift increases non-linearly. There- 
fore, out-of-plane excursion is particularly harmful. As 
long as the magnetization vector is not deflected too far 
out of the magnet's plane, and does not stray into the 
ambit of another stable state, it always returns to the 
original state and recovers. This self-correcting behavior 
lends itself to image recovery. 

Fig. 2 (a) shows the scheme for encoding pixel shades 
in a black/gray /white image. The magnetization's orien- 
tation encodes three different shades as follows: ^ = 0° 
(black), (/) = 90° (gray), (j) = 270° (gray) and (/) = 180° 
(white). The shades can also be assigned numerical val- 
ues: black = 0, gray = 0.5 and white = 1. 

A 512x512 pixel image is encoded by the above scheme 
and is shown in Fig. 2 (b). Next, the numerical value 
of each pixel is changed randomly to simulate the effect 
of noise. We restrict the random out-of-plane deflection 



of the magnetization vector to ±1°, which then re- 
stricts the in-plane deflection Acj) to ±40° since that is 
the maximum in-plane deflection that can be corrected 
when |A^| < 1°. The choice of ±1° out-of-plane deflec- 
tion is dictated by the fact that this allows a reasonably 
large azimuthal deflection. This visually distorts the im- 
age by a large degree in Fig. 2(d). In accordance with 
this choice, a pixel with intensity value = is randomly 
assigned a value between and 0.222. If the intensity 
value = 0.5, it is changed to something between 0.278 
and 0.722. Similarly, for pixels with intensity value = 
1, the value is changed to something between 0.778 and 
1. These distortions restrict the azimuthal deflections to 
±40°, i.e. |A0| < 40°. 

Fig. 2 (c) shows the corrupted image at t = ns for a 
critical region (the bird's eye). By converting each pixel 
of the corrupted image to its equivalent ^- value and then 
solving the coupled equations for 6 — (j) dynamics^, we 
can determine the final state (black, gray or white) of 
the pixel. Figs. 2 (d), (e) and (f) illustrate the image 
recovery process. At t = 1 ns, the noise in the image has 
been greatly reduced, with most pixels settling back into 
their original states (Fig. 2 (e)). Steady state is achieved 
in 2 ns (Fig. 2 (f)) since the images at 2 ns and at 3 
ns (the latter not shown here) are identical. The final 
steady state image is identical with the original image 
pixel by pixel, showing 100% recovery. 

One could extend this scheme to pattern recognition as 
well. This would be achieved by integrating a magnetic 
tunnel junction (MTJ) vertically underneath each mag- 
net. The magnetization orientations of the lower hard 
magnetic layers store the image to be compared against. 
The input image is written in the upper soft magnetic lay- 
ers (Ni). If the input matches the stored image pixel by 
pixel, then the magnetizations of the soft and hard mag- 
nets will be all parallel. This will result in the maximum 
tunnel current flowing through each MTJ. By setting an 
appropriate current threshold (say X% of the maximum) , 
we can determine if the two images match with probabil- 
ity of X% and thus recognize the input image. 
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